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32:   Lines  12  -  14  should  be  replaced  by  the  following  sentence 
"It  was  pointed  out  in  the  previous  section  that  when  the 
applied  force  varies  sinusoidally  in  time,  the  fields  in 
the  plasma  also  vary  sinusoidally  in  time  with  constant 
amplitude  except  at  the  point  x  =  x  (w  )  where  the  fields 
grow  in  time  . ,T 

36:   In  line  13,  the  number  ''0.07''  should  be  replaced  by  "7." 


After  submission  of  the  report  for  printing,  the  following  referei 
concerning  Alfven  wave  heating  was  referred  to  the  author: 
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Abstract 

The  relationship  between  energy  absorption  and  the  continu- 
ous frequency  spectrum  of  the  linearized  equations  of  ideal  mag- 
netohydrodynamics  Is  investigated.   We  limit  our  considerations 
to  incompressible  fluid  perturbations  for  which  the  continuum 
stems  from  resonant  surfaces  where  the  oscillation  frequency 
equals  the  local  frequency  of  an  Alfve"n  wave.   Details  of  this 
absorption  process  are  illustrated  by  obtaining  the  response 
of  the  diffuse  sheet  pinch  to  an  external  current  source. 
Expressions  are  derived  which  give  the  rate  at  which  energy  is 
transferred  to  the  plasma  and  show  the  spatial  distribution 
of  the  absorbed  energy.   The  results  indicate  that  the  absorp- 
tion is  enhanced  if  the  plasma  profile  has  regions  of  large 
spatial  gradients.   This  enhanced  absorption  is  a  consequence 
of  the  presence  of  surface  waves. 


1^  Introduction 

Conservative,  dynamical  systems  under  the  linear  approxi- 
mation have  the  property  that  they  can  display  dissipative 
effects  if  there  exists  a  sufficiently  large  number  of  degrees 
of  freedom  (Sturrock  1961;  Grad  1969).   In  simplest  terms 
apparent  dissipation  arises  when  the  system  has  a  continuous 
frequency  spectrum.   The  normal  modes  which  make  up  the  con- 
tinuum can  phase-mix  in  time,  leading  to  decay  of  macroscopic 
variables.   Plasma  models,  free  of  any  apparent  loss  mechanism, 
are  often  encountered  which  are  characterized  by  a  continuum 
and  are  consequently  dissipative.   Among  these  models  we 
mention  a  thermal  plasma  described  by  the  linearized  Vlasov 
equation  (Landau  1952;  Weitzner  1963),  a  cold  non-uniform 
plasma  (Barston  1964)  and  certain  configurations,  such  as  the 
screw  pinch  (Grad  1973;  Goedbloed  and  Sakanaka  197^),  the 
theta  pinch  (Weitzner  1973;  Grossmann  and  Tataronis  1973)  and 
the  sheet  pinch  (Boris  1968;  Tataronis  and  Grossmann  1972; 
Uberoi  1972;  Tataronis  and  Grossmann  1973),  described  by  the 
linearized  equations  of  ideal  magnetohydrodynamics .   (In  what 
follows  the  initials  MHD  will  designate  the  term  magnetohydro- 
dynamics) . 

A  second  implication  of  a  continuum  is  that  the  system 
can  irreversibly  abosrb  energy  when  under  the  influence  of  an 
externally  applied  force.   The  absorption  is  a  consequence  of 
a  resonance  which  occurs  when  the  spatial  distribution  and  the 
time  behavior  of  the  force  matches  a  mode  in  the  continuum. 
A  characteristic  of  this  absorption  is  that  it  occurs  locally 


in  the  space  in  which  the  continuum  is  defined.   This  relation- 
ship between  a  continuous  spectrum  and  energy  absorption  has 
been  recognized  in  the  past  (Grad  1969;  Baldwin  and  Ignat  1969; 
Tataronis  and  Grossmann  1972;  Hasegawa  1973),  and  calculations 
have  been  made  to  determine  the  effectiveness  of  this  method 
of  energy  transfer  as  a  means  of  heating  plasmas  (Baldwin  and 
Ignat  1969;  Brownell  1972;  Tataronis  1973;  Grossmann,  Kaufmann 
and  Neuhauser  1973)- 

In  what  follows  we  shall  be  concerned  with  this  absorption 
mechanism  under  the  approximation  of  ideal  MHD.    We  shall 
limit  our  considerations  to  a  particular  plasma  configuration, 
namely  the  diffuse  sheet  pinch,  for  which  it  has  already  been 
shown  that  the  linearized  MHD  equations  have  a  continuous  spec- 
trum (Boris  1968;  Tataronis  and  Grossmann  1972;  Uberoi  1972; 
Tataronis  and  Grossmann  1973)-   Under  the  approximation  that 
the  fluid  motion  is  incompressible,  we  obtain  the  response  of 
the  plasma  to  an  external  sheet  current  located  in  a  surroundinf 
vacuum  region  and  determine  from  the  response  the  energy 
transferred  to  the  plasma.   The  force  which  drives  the  plasma 
is  derived  from  the  pressure  exerted  on  the  plasma  surface 
by  the  magnetic  field  induced  in  the  vacuum  by  the  sheet  cur- 
rent and  the  plasma  current.   As  we  shall  demonstrate,  this 
scheme  leads  to  efficient  coupling  of  energy  to  the  interior 
of  the  plasma  region  even  though  the  plasma  motion  is 


The  results  of  this  paper  were  first  reported  at  the  October 
1973  meeting  of  the  Plasma  Physics  Division  of  the  American 
Physical  Society,  Philadelphia,  Pa. (Abstract  in  Bull.  Am.  Phys . 
Soc.  18,  1291  C1973)). 


incompressible.   The  absorption  rate  is  measured  in  terms  of 
the  effective  impedance. that  the  plasma  imposes  on  the  sheet 
current.   It  will  be  shown  that  as  a  consequence  of  the 
continuum,  the  impedance  has  a  resistive  component  which,  for 
parameters  and  dimensions  typical  of  high  beta,  high  density 
plasmas,  can  have  values  exceeding  one  Ohm  at  frequencies  on 
the  order  of  2  MHz.   This  resistance  is  comparable  to  what 
has  been  found  for  excitation  of  the  lower  hybrid  resonance 
by  a  method  similar  to  the  one  assumed  here  (Puri  1973) >  with 
the  advantage  that  the  absorption  described  below  occurs  at  a 
relatively  low  frequency. 


2.   The  Plasma  Model  and  Basic  Equations 

The  specific  plasma  model  we  study  is  the  planar  sheet 
pinch  with  the  equilibrium  magnetic  field  J3  C=  B  z)  parallel 
to  the  z  axis  of  an  (x,y,z)  Cartesian  coordinate  system  as 
illustrated  in  Fig.  1.   The  equilibrium  pressure  p(x)  and 
mass  density  p(x),  as  well  as  B^Cx),  depend  only  on  x  and  are 
symmetric  with  respect  to  x  =  0.   The  plasma  is  confined  to 
the  region  -b <  x  < b .   For  |x|  > b ,  there  is  a  vacuum  where 
p  =0,  ^  is  independent  of  x  and  an  externally  supported  sheet 
current  sheet  £      is  present  which  is  directed  parallel  to  the 
y  axis  and  depends  sinusoidally  on  z, 

£s(x,y,z,t)  =  JsCt)[SCx-c)  -  5tx+c)]cos  kz  y         CD 

In  this  expression  k  is  a  real  wave  number,  6Cx)  Is  the  Dirac 
delta  function  £  has  the  dimensions  of  current  per  unit  area 
and  J  (t)  is  expressed  in  terms  of  current  per  unit  length, 
length  being  measured  in  the  z  direction.   The  sheet  current 
is  the  source  of  the  external  force  which  acts  on  the  plasma 
surface. 

We  descrihe  the  plasma  dynamics  by  means  of  the  linearized 
equations  of  ideal  MHD  with  resistivity  and  viscosity  neglected 
and  incompressibility  assumed.  In  Maxwell's  equations,  the 
displacement  current  is  neglected.  Since  the  assumed  equili- 
brium magnetic  field  is  free  of  curvature,  we  have  (Jg  •  V)£  =  0, 
implying  that  the  MHD  equations  can  be  reduced  to  the  following 
form: 


vi£  (2) 

at' 


P^-^->V-1f 


.g-0  (3) 


CO 


^  -  -V  x  ;  C5) 


*1  =  £  V  x  8i  C6) 


where  y  is  the  vacuum  permeabily,  ^,  jgp  and  £  are  perturbations 

of  the  velocity,  magnetic  and  electric  fields,  and  p   represents 

*         —1 
a  perturbation  of  the  total  pressure,  i.e.  p   =p1  +  y  £*$]_• 

In  the  vacuum  the  fluctuating  electric,  6£v,  and  magnetic,  6£y, 

fields  are  determined  by, 


6Bv  =  Vc|> 
3( 


rr  -  "v  x  6Sv 


C7) 


(8) 


where  the  magnetic  potential  <J>  satisfies, 

V2(j>  =  0.  (9) 

At  x  =  b,  the  plasma  and  vacuum  fields  are  connected  by  the 
usual  continuity  conditions  of  ideal  MHD,  namely  6By  is  parallel 
to  the  fluctuating  plasma-vacuum  interface  and  the  total  pres- 
sure p  + B2/u  is  continuous  across  the  interface.   After 


linearization,  these  two  conditions  imply  the  following  equali- 
ties at  x  =  b, 

8  x  6Sv  =  CQ  "«^v  Cl0) 


jl  £v  '  6£v  =  P*  CXI] 


where  B   (  =  B   2)  is  the  equilibrium  magnetic  field  in  the 
vacuum  region  and  (i  is  a  unit  vector  normal  to  the  plasma- 
vacuum  interface  of  the  equilibrium  and  directed  into  the 
vacuum.   Finally,  the  following  jump  conditions  across  the 
sheet  current  are  to  be  satisfied, 

n,   x    C6^v]    =   y^s  C12] 

e  -  isj^j  =  o  ci3i 

where  the  brackets  designate  the  jump  in  the  enclosed  quantity. 
The  symmetry  of  the  sheet  current  implies  that  the  excited 
mode  is  characterized  as  follows, 


v,  =  (vx,  0,  vz),  (14) 


1  =  CO,  Ey,  0),  (15) 

El  =  CBx'  °-  3z}>  (16) 

h   =  (o,  V  0)j  Cl7) 


6£v  =  (6Bx,  0,  6Bz),  (18) 

6£v  =  (0,  6Ey,  0)  (19) 

where  in  each  expression  the  parenthesis  on  the  right  enclose 
the  x,  y,  and  z  Cartesian  components  of  the  vector  on  the  left 
Furthermore,  it  can  be  readily  established  that  the  Induced 
fields  possess  a  spatial  symmetry  with  respect  to  the  plane 
x  =  0,  namely, 


vx(-x,z,t)  =  -vxCx,z,t),  vz(-x,z,t)  =  vz(x,z,t),      (20) 

Ey(-x,z,t)  =  -Ey(x,z,t),  (21) 

Bx(-x,z,t)  =  -Bx(x,z,t),  Bz(-x,z,t)  =  Bz(x,z,t),       (22) 

Jy(-x,z,t)  =  -Jy(x,z,t),  (23) 

p*(-x,z,t)  =  p*(x,z,t).  (2H) 


SB   and  6E  have  the  same  symmetry  as  B,  and  E  respectively. 

Finally,  since  the  imposed  surface  current,  eq.  (1),  depends 

on  z  in  the  form  cos  kz,  it  follows  that  v  ,  B  ,  E  ,  J  ,  p*, 

SB   and  6E   also  have  the  form  f(x,t)  cos  kz,  while  v  ,  B   and 
z       y  z    x 

SB  depend  on  z  as  g  (x,t)  sin  kz.   In  what  follows  the  depen- 
dence of  the  fields  on  z  as  either  coz  kz  or  sin  kz  will  not  be 
explicitly  shown  except  in  cases  where  ambiguity  would  result 
by  leaving  the  z  dependence  implicit. 


In  the  sheet  pinch  geometry,  the  plasma  fields  and  energy 
can  be  expressed  explicitly  in  terms  of  v  .   In  particular, 
by  combining  eqs .  (2)  and  (3)5  one  finds, 

1  3vx 
Vz  =  ~t  IT  (Xjt)  sln  kz>  (25) 


^  =  -■-l£fiCxJ.4T  +  £  k2B2(x)]^Cx,t)cos  kz,  (26) 


where  v  satisfies, 

gk(  -\  +  i-  k2B2)^]  -k2(p-^2  +  h   t2B2)Y-  =  0.       (27) 
3x    3t2    y        3x         9t2    y         X 


The  remaining  fields  are  obtained  from  v  and  v  with  the  aide 
of  eqs.  (4)-(6).   As  a  consequence  of  the  spatial  symmetry  in- 
dicated by  eqs.  (20)-(2il),  it  suffices  to  solve  eq.  (27)  in 
the  region  x  >  0  subject  to  boundary  conditions  at  x  =  0 , 


vx(0,z,t)  =  0,  (28) 

and  at  x  =  b , 

v  (b,z,t)  =  jL  6E(b,z,t),  (291 

x  v   y 

p*(b,z,t)  =  ^  Bv6Bz(b,z,t),  (30) 


where  <5E  (x,z,t)  and  6Bz(x,z,t)  satisfy  eqs.  (7)-(9),  are 
consistent  with  eq.  (12)  and  (13) >  and  vanish  as  x  ■*  od  . 

The  equation  which  expresses  the  rate  of  change  of  the 
plasma  energy  in  terms  of  v  is  derived  by  introducing  ^  (x,z,t), 
the  displacement  of  a  fluid  element  from  its  equilibrium  position 


3 


t 

£  =  f  «•*  (3D 

0 


By  integrating  eq.  (2)  with  respect  to  t,  with  the  initial 
conditions  set  equal  to  zero,  scalar  multiplying  by  y_  and 
the  integrating  over  some  plasma  volume  V  bounded  by  a  sur- 
face E  ,  one  finds  the  following  statement  of  energy  conserva- 
tion in  terms  of  ^  and  %, 
2 

l  d&¥   TT  "  ?  X,  '  ^  •  v>2^  -  -  f  «*  ■  *P-  (32) 

V  £ 

where  use  was  made  of  the  incompressibility  condition  V  •  y/  =  0 
For  the  volume  V,  we  choose  the  region  bounded  on  the  x  axis 
by  x  =  ±b ,  on  the  y  axis  by  y  =  0  and  y  =  L  ,  where  L   is  an 


arbitrarily  chosen  dimension,  and  on  the  z  axis  by  z  =  0  and 
z  =  2ir/k  .   By  taking  into  consideration  that  v  and  p*  are 
antisymmetric  in  x  and  depend  on  z  in  the  form  cos  kz  while 


integration  and  the  integration  over  y  and  z  in  the  left  hand 
side  of  eq.  (32)  can  be  readily  accomplished,  yielding, 

dT  =  -"IT2"  vx(b,t)p*(b,t),  (33) 

where  W,  the  plasma  energy  ,  is  given  by, 

W  =  -^Z-   f  dx[|p(x)^2(x,t)  +  Jp  k2B2(x)£2(x,t)]  .        (3<0 

Q 

2  2     2  2 

^   signifies  v   +  v  ,  with  a  similar  definition  for  ^ 


Equation  (33)  states  that  the  rate  at  which  the  plasma 
energy  increases  in  the  enclosed  volume  is  equal  to  the  energy 
flux  density,  v  (b  ,t )p*(b ,t ) ,  averaged  over  y  and  z  at  x  =  ±b , 
multiplied  by  the  surface  area  2ttL  /k.   It  follows  from  the 
boundary  conditions  given  by  eqs .  (29)  and  (30)  that  the  energy 
flux  density  is  equal  to  the  x  component  of  the  vacuum  Poynting 
vector  S   e   <5E   x  &%v/v  at  the  plasma  boundary,  which  in 

turn  can  be  related  by  means  of  Maxwell's  equations  to  the  rate 
at  which  the  energy  of  the  perturbed  vacuum  magnetic  field  W 
changes  in  time  and  to  the  power  delivered  by  the  sheet 
current,  -(2tiL  /k)  6E  (c,t)  Jg(t).   Explicitly,  eq.  (.33)  is 
equivalent  to  the  following  expressions, 


,w  2ttL      6E(b,t)<SB    (b,t) 

ill  =    _   L  — I Z (l^) 

dt  k  u  LJ5; 


2irL 


=   "  ~dT  "  ~k       6Ey(c'^    JsCt),  (36) 

where, 

2ttL       °°r  6B2 

Wv    =    V-     J    dX    ^    (Xjt)  (37) 


Equations  (.25) -(30)  and  (33)-(37)  form  the  basis  of  the 
remaining  sections  of  the  paper.   In  section  3  we  solve  eq. 
(27)  for  v   and  use  this  solution  to  obtain  expressions  for 
the  plasma  and  vacuum  fields  in  terms  of  the  surface  current. 
Section  4  is  devoted  to  an  examination  of  the  temporal  behav- 
ior  and  spatial  distribution  of  the  energy  absorbed  by  the 
plasma,  and  finally  in  section  5  we  examine  the  absorption  rate 
for  a  sp0"^'"  Dlasma  profile. 


In 


3.   Solution  of  MHD  Equations 

In  this  section  we  make  use  of  the  Laplace  transform, 
defined  by, 

vx(x,t)  =  j  |^-e1(l)t  tfx(x,w)  ,    Im(w)<  0  ,         (38) 


vx(x,w)  =   J  dte"lwt  vx(x,t)  , 


(39) 


to  solve  the  linearized  MHD  equations  for  the  plasma  and  vacuum 
fields.   In  eq.  (38),  *£  is  a  contour  of  integration  below  the 
singularities  of  the  integrand  in  the  lower  half  complex  w 
plane.   The  differential  equation  which  v   satisfies  is  obtained 
by  transforming  eq.  (27)  with  eq.  (39),  yielding, 

^r[A(x,03>^]  -  k^A(x,co)v   =  0  ,    Im(u>)  <  0  ,  (40) 

oX       '        a  X  X 


where , 


A(x,u>)  =  -p(x)u)2  +  k2B2(x)/y  ,  (4l) 


and  initial  conditions  have  been  set  equal  to  zero.  From 
eqs.  (25)  and  (26),  the  Laplace  transforms  of  vz(x,t)  and 
p*(x,t)  are, 


,  9v 
Vx»u*  =  -i^Cx.o,)  ,  (42) 


3V 
iwk^   ax 


11 


The  vacuum  fields,  defined  by  eqs .  (7)-(9),  can  be  explicitly 

+  kx 
expressed  in  terms  of  the  exponential  function  e    .   In 

particular,  the  vacuum  electric  field  which  satisfies  the 

jump  conditions  at  the  sheet  current,  and  vanishes  as  x  •>  oo  , 

is,  after  performing  a  Laplace  transform, 


yJ  (w) 
■iu — n e  ""^  sinh  kx  +  D(oj)e  ^  ,   b<  x<  c  (44) 


6Ew(x,tai) 

-iu — r e"^  sinh  kc  +  D(u)e~*-A  ,   x>  c    (45) 


,  »JaM    .-kx  _ 


where  D(u)  is  a  function  of  u  determined  by  the  boundary 
conditions  at  x  =  b .   The  electric  field  has  been  written  as 
the  sum  of  two  terms:   the  first  term  represents  the  electric 
field  induced  in  the  vacuum  directly  by  the  sheet  current,  while 
the  second  term  Is  the  component  of  the  electric  field  resulting 
from  the  plasma  motion.   If  the  plasma  is  removed,  then  the 
second  term  vanishes.   The  boundary  conditions  at  x  =  b ,  eqs. 
(10)  and  (11),  can  be  expressed  in  terms  of  v  and  E  .   After 
a  Laplace  transformation  is  performed,  they  read, 


v(b,w)  =  =±-   SE(b,w)  ,  (46) 

x        By   y 


9v  k2B    36E 


(47) 


At  x  =  0  we  have 


v  C0,w)  =  0  .  (48) 


12 


Equations  (40 ) — C^8 )  form  a  self-consistent  system  which  deter- 
mines v  and  E  . 
x      y 

To  solve  the  system  we  introduce  a  function  \\>    (x,cd)  which 
satisfies  eq.  (40),  is  free  of  singularities  on  the  real  x 
axis  in  the  range  0  _<  x  _<  b  for  w  on  the  Laplace  contour,  and 
satisfies  the  following  conditions  at  x  =  0 , 

K0,w)  =  0  ,  (49) 


§f(0,u>)  =  k  .  (50) 


v  (x,w)  must  then  be  proportional  to  i> : 


vxCx,w)  =  C(a>)Kx,uO  (51) 


The  constants  C(w)  and  D(w)  are  obtained  by  substituting  eqs . 
C44)  and  (51)  into  eqs.  (46)  and  (47)  yielding  for  lm(w)<0, 

uJs(w)   k2B2       e-k(c-b) 
C(o))  =  -io)   Bv    ^7A(b7^  Ll+kq(uO]<K(b,w)  (52) 

yJ-^w^   °,inh  kb  -  kq(o))  cosh  kb   -k(c-b)       ,„% 
1  +  kq(u>)        e         '     °J; 


where  iJj*(x,u))  signifies  d^(x,u)/dx  and, 

k2B2 

v    J(^4   .  (51,) 


q(u)  E  pACb,o))  #urS~ 


It  should  be  noted  that  q  and  v  are  independent  of  the  normal- 
ization imposed  on  ty. 
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i>   can  be  expressed  in  terms  of  known  functions  only  in 
certain  cases.   One  example  is  the  case  where  A  is  independent 
of  x,  for  which  eq.  (40)  has  the  solution, 

4>   =  sinh  kx,  (55) 

A  second  example  is  the  piecewise  smooth  profile  shown  in 
Fig.  2.   In  this  case  the  plasma  equilibrium  is  uniform  in 
the  region  |x|  <a,  while  in  the  sheath  a<  [x[  <b,  p(x)  and  B  (x) 
vary  linearly  with  x,  implying, 

2    2 

2  k  K 

a)   p     +   £— ,        -a   <  x    <  a   , 

Mp  M  -       - 

(56) 

w2pb_^xl+   ki[B2  2    _   B2)Mib] 

Mp  b-a         y  v  v  py    b-a   J 


<  x    <  b 


and, 


sinh  kx,  0    <  x    <   a, 


p(U)I[kCxo-x)]-Ko[k(x-x)]  (57) 

Slnh  ka     p^)to[kCxo-a)]-Kork(xo-a)]    -    a<x<b 


where  I  (z)  and  K  (z)  are  modified  Bessel  functions  of  order 
zero  as  defined  by  Watson  (Watson  1958),  x  Cu>)  is  defined  by 


A(x  ,oj)  =  0,  and. 


KQ[k(xo-a)]  -  K1[k(xQ-a)]tanh  ka 
p(a))  =  I°[k(xo-a)J  +  I1Lk(xQ-a)]tanh  ka   ' 


At    x    =    a,    we    imposed   the  constraint  that    ijj/iK    be    continous 


Ik 


More  generally,  if  A(x,u))  posesses  a  power  series  expansion 
in  some  domain  about  x  (u)  which  includes  0  <  x  <  b, 


o 


A(x,w)  =JZ1   an(x-xjn  ,  (59) 


n=l   n 


then  for  a,  ?   0,  i.e.  for  A'(x  ,w)  /  0,  it  can  be  shown  that 
about  x  ,  two  independent  solutions  of  eq.  (40)  are, 


f-L  =  g(x,ul  ,  (60) 


f2  =  g(x,a))  9m    Cxq-x)  +  (x-xo)h(x,w),  (61) 


with  which  we  form, 


*  =  p(ai)g(x,u)  +  EgCx,u)  V*j    (xQ-x)  +  (x-xo)h(x,a))],        (62) 


p(cj)  is  independent  of  x,  and  g(x,co)  and  h(x,w)  are  functions 
analytic  in  some  domain  about  x  (w),  the  extent  of  the  domain 
being  determined  by  the  closest  zero  of  A(x,u)  to  x  .   By 
analytic  continuation,  this  domain  can  be  extended  to  include 
x  =  0  where  the  boundary  conditions  (49)  and  (50)  are  imposed. 
This  uniquely  determines  g,  h  and  p. 

Both  eqs.  (57)  and  (62)  show  that  when  A  depends  on  x,  i> 
possesses  a  logarithmic  singularity  at  the  point  where  x  (to)  =  x. 
Since  the  origin  and  significance  of  this  singularity  are 
considered  elsewhere  (Tataronis  and  Grossmann  1973;  Grossmann 
and  Tataronis  1973),  we  shall  limit  our  comments  here  concerning 
it  by  mentioning  that  the  singularity  can  be  related  to  the 


15 


existence  of  a  continuous  spectrum  in  the  frequency  range  where 
A(x,oo)  vanishes  at  some  point  x  (to)  in  the  plasma  region. 
From  eq.  (41),  the  continuum  has  the  following  bounds, 


2         2         2 
min[to  (x)]  <  to   <  max[u  (x)]  ,  (63) 


2        2  2 
where  u  (x)  =  k  B  (x)/up(x)  is  the  square  of  the  local  fre- 
quency of  an  Alfven  wave  propagating  with  wave  number  k.   If 
to  satisfies  eq.  (63),  then  there  exists  a  point  x  (to)  in  the 


plasma,  defined  by 


w2  =  <?(x  )  ,  (64) 


ularity.   By  differentiating  v  with  respect  to  x  to  form  v 
according  to  eq.  (42),  one  finds  that  v   is  also  singular  about 

x  > 
o' 


v  ^  1^  (65) 

z    x-x  (to)  ' 


but  the  singularity  is  stronger  than  that  of  vx  in  the  sense 
that  when  oj  is  in  the  continuum,  v  is  not  square  Integrable 


integrable  which  implies  the  absorption  process  that  is  the 
subject  of  the  present  paper. 

The.  presence  of  a  logarithmic  factor  renders  ty  a  multi- 
valued function.  Branch  cuts  must  then  be  introduced  in  the 
complex  co  plane  such  that  with  to  on  the  principal  branch, 
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ty  is  free  of  singularities  on  the  real  x  axis  in  the  range 

-b <  x  <b.   For  the  piecewise  smooth  profile  in  Fig.  2,  4,  is 

given  by  eq.  (57),  and  the  associated  branch  points  and  branch 

cuts  are  shown  in  Fig.  3.   On  the  positive  real  axis  of  the 

w  plane,  ^  is  real  for  w  >  w(x).   Due  to  the  logarithmic 

factor  in  K  ,  a  branch  point  is  present  at  w  =-u_(x),  with 
o  a 

the  corresponding  branch  cut  defined  by, 


a)  -  03  (x)  =  I  a)  -  u)  (x)|e   ,    TT  >  6  >  — rr  •  (66) 


It    can  be   shown   that  for  xq    >  a,    <*q.    (66)    implies 


i8 
x    (w)    -   x   =    |  xo(o>)    -  x|  e       ,        v    >    B    >    ~n ,  C67a) 

[p(co)+i€7r]I    [k(x    -x)]-K    [k|x   -x|] 

>  =   sinh  ka  * r    w  r°, ?.  v   r r°, °  (67b) 

p(o))I  [k(x  -a)  j-K  Lk(x  -a)] 


where  to  i-s  real,  and  €  p  +1,  -l  or  0  for  6  =  +rr,  -it  or  0. 
Therefore,  along  the  branch  cut,  i|>  is  complex  with  the  sign 
of  its  imaginary  part  dependent  on  whether  w  lies  on  the 
upper  (6=tt)  or  lower  (0=  -v)    edge  of  the  cut.   The  branch  cut 
terminates  at  the  second  branch  point  uo  =  00  .   For  0  <  qj  <  u  , 
ty   is  again  real.   The  behavior  of  ^  on  the  negative  w  axis  is 
given  by  the  symmetry  relation, 


(x,  -co)  =  /(x,u)),  (68) 


where  y     designates  the  complex  conjugate  of  4>.   As  we  will 
show  in  the  next  section,  the  imaginary  part  of  ty   plays  an 
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important  role  in  demonstrating  that  the  plasma  absorbs 

energy  when  the  applied  frequency  is  in  the  continuum. 

p 
If  to „(x)  is  an  analytic  function,  one  can  construct  the 

branch  cuts  of  ^  by  using  either  eq.  (62)  or  more  conve- 
niently the  following  representation  of  4> } 

*  =  Q(u)[p(w)*1Cx,u)  +  <|/2(x,w)],  (69) 

where  Q  and  p  are  functions  of  co  determined  by  the  boundary 
conditions  at  x  =  Q,  and  ij;,  and  ty-   are  two  linearly  independent 
solutions  of  eq.  (40)  which  we  write  in  the  following  form, 

*1  =  g(x,a3),  (70) 


dx' 


i>2   =  g(x,u>)   j  - 


ACx>)g"cx;co) 


g(x,to)  is  analytic  about  x   and  is  assumed  to  have  the  norm- 
alization g(x  ,to)  =  1.   For  a  particular  value  of  x,  x  =  x  , 
the  branch  cuts  of  4>    (x  ,to)  coincide  with  those  ranges  of  go 
in  which  A(x,co)  vanishes  somewhere   for  0  <x<x  .   The  branch 

points  are  located  at  oj  =  ±oj  (x  )  and  at  those  values  of  to  (x) 

a  P  a 

p 
at  which  dco  /dx  =  0.   To  establish  the  existence  of  a  branch 
a 

2 
point  at  co  (x  ),  we  expand  w  (x)  about  x  , 

.  2 
</(x)  =  .2a(xp)  +  ^(Xp)(x-xp)  +  ...  (72) 


to  derive  the  following  expansion  of  x  (to)  about  to  (x  ) 
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2   2,   n 

10  -w  (x  ) 

x  Coo)  ^  xr  +  a  P  -  +  ...  (73) 

o       p    (duT/dx), 

Assuming  (dto2/dx)Y=Y  >  0,  eq.  (73),  together  with, 
a    *  xp 

oj  -  to  (x  )  ■  |(rt  -  to  (x  )|ei9,    tt  >  9  >  -it        (74) 
3.  p  a.  p  —   —  ■ 


implies   the    following  behavior   of   \\>2 (x    ,to)    on   the   real   oj   axis 

near   to    (x    )  , 
a     p 

x 

*P(x    .u>)    =    gCxn)oj)P        ^p r  +   2^ '    (75) 

2      P  P  I      ACxiuJg^Cxiu)       p(x    UdoZ/dx) 

d  p       a.       A  — A 

where  P  designates  that  the  principal  value  of  the  integral  is 

to  be  taken  when  x   lies  on  the  real  x*  axis,  and, 
o 


(76) 


Equation  (75)  implies  the  existence  of  a  branch  cut  on  the 

real  axis  for  to  >  to  (x  )  .   Since  the  boundary  conditions  which 
a  p 

iJj  must  satisfy  at  x  =  0  are  real,  it  is  readily  established 

with  eq.  (75)  that  Kx  ,to)  is  real  for  to  >  <*>a(x  )  while  complex 

2 
along  the  branch  cut.  A  similar  conclusion  follows  for  (d(oa/dx)x  « 

2  P 

Let  us  now  assume  that  dto  /dx  =  0  at  x  =  x   and  that  about 

x  ,  ^2(x)    can  be  expanded  in  the  following  manner, 


w2(x)  =  to2  +  Bv(x-xv)2  +  ...  C77) 
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Equation    (77)    implies   that    for  to  %   u>v ,   A(x,co)    possesses   two 
zeros,    x    i(iii]    and  x   2Cto)  ,   which,  have   the   following  expansions 


2      2 
03    -o)v    ,  /0 

xvlCoi)  *  xv  +  C— g^)1/2     +  ... 

2     2  C78) 

XV2<»    *   xv    "    C~T^)  V        +    ••' 

2 
Assuming  6   >  0,  corresponding  to  a  local  minimum  of  w  (x), 
v  a 


eq.  (78)  combined  with, 


to  -  tov  =  |u  -  o)v|e1(e"27r)  ,   0  <  6  <  2tt         (79) 


implies   for   xv   <    a  <   x      and  w  real, 


<K(x.co)    =    g(x_,u)P    f    dX'        ,  ^^T1  DfxY 

2      P  P  i    A(x;o))g2Cx;a))  ^     PUp}Bv 


where , 


6„      1/2 


C-^2) 


(80) 


(81) 


A  similar  expression  can  be  derived  for  f,  if  6   <  0,  corres- 

2 
ponding  to  a  local  maximum  of  to  (x) .   The  implication  of  eq. 

(80)  is  that  4>   possesses  a  branch  point  at  to  ,  with  the  branch 

cut  extending  to  the  region  to  >  u  .   Figure  H    indicates  the 

branch  cuts  which  are  implied  by  the  above  considerations  when 
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to  (x)  depends  on  x  in  the  manner  shown  in  Fig.  4a.   It  can 

2 
be  shown  that  when  <o  (x)  Is  analytic,  eq.  (68)  is  satisfied 

for  real  co. 

We  now  use  the  explicit  forms  of  i|>  derived  above  to  carry 
out  the  Inverse  Laplace  transform  of  v  (x,to)  to  obtain  v  (x,t). 
In  eq.  (38)  the  integration  on  the  Laplace  contour  o£   can  be 
performed  be  exploiting  the  singularities  of  the  integrand  in 
the  complex  u  plane.  The  procedure  we  follow  is  standard  in 
the  theory  of  Laplace  transforms,  namely,  we  deform  the  contour 
/from  its  position  in  the  lower  half  <o  plane  below  the  singu- 
larities of  v  (x,t)  to  a  contour^*  which  enciricles  the  singu- 
larities and  extends  into  the  upper  half  plane.   Even  if  the 
integration  cannot  be  explicitly  carried  out  with  the  aide  of 
this  new  contour,  one  can  obtain  with*£  a  time  asymptotic 
expansion  of  the  integral  be  expanding  the  integrand  about 
the  singularities  and  then  performing  the  integration  for  large 
t.   We  follow  this  procedure  for  the  case  where  the  sheet 
current  varies  sinusoidally  in  time, 


J  (t)  =  J  sin  u  t,   J  (co)  =  J  -x-By   •  C82) 

s        o       o      s        o  *_^ 

o 

After  combining  eqs.  (38),  (51),  (52)  and  (82)  and  deforming 
jC  to  oC    in  the  upper  half  plane,  one  obtains  for  ioq  /   <oa(x), 


v  (x,t)  =  JoRv(x,coo)  sin  uQt  +  JoXv(x,coo)  cos  ujb  +  vxt(x,t)  .  (83) 


R  and  X   are  real  functions  defined  by 


21 


R      +    i    LfE    Z       ,  (84) 

where , 

v  .  k2B2  -k(c-b)    .  ,         x 

Zv(x,o))    =    J-  =    -  -B^uA(b,o))    I+l^To)!  HG)       '  C85) 

s 


=    j      §|  eiwt    Z 


xtv    '    '  2tt  -vv"'-'    ~s 


As  illustrated  in  Fig.  5,^  extends  into  the  upper  half  plane 
around  the  singularities  of  Z   .   Z   is  characterized  by  two 
sets  of  singularities.   The  first  set  consists  of  the  branch 
points  and  corresponding  branch  cuts  of  ty   on  the  real  axis. 
In  Fig.  5  the  branch  cuts  have  been  analytically  continued 
into  the  upper  half  plane.   By  performing  this  analytic 
continuation,  one  obtains  the  second  set  of  singularities, 


/LCw  )  =   1  +  kq(co  )  =  0  .  (87) 

s  o 


It  is  pointed  out  that  under  the  incompressible  approximation 
A-  does  not  possess  zeros  on  the  principal  sheet  of  the  w  plane. 
Zeros  are  found  only  after  an  analytic  continuation  has  been 
performed  across  the  logarithic  branch  cut.   The  significance 
of  the  roots  of  eq.  (87)  are  discussed  in  considerable  detail 
elsewhere  (Tataronis  and  Grossmann  1973;  Grossmann  and  Tataronis 
1973).   Therefore,  we  shall  limit  our  comments  here  concerning 
them  by  indicating  that  w   is  in  general  complex  and  corresponds 
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to  the  frequency  of  a  surface  wave  which  decays  expontially 
in  time  due  to  phase  mixing.   This  decay  is  analagous  to  the 
decay  of  space  charge  waves  in  a  Vlasov  plasma  by  the  mechanism 
of  "Landau  damping"  . 

As  t  -*■  oo  ,  one  finds 

COS  to  (x)t  COS  to   t 

vxt(x,t)  %  aCx)      (;)t  +CYvCx) 


(a>  t) 
v  av 


(88) 


+  S(x)e-U,sit  cos  tot 


where  co   and  w  .  are  the  real  and  imaginary  parts  of  to  ,  and 
a,  Yvand  S  are  functions  of  x  obtained  by  expanding  Z  (x,a>) 
about  each  singularity  and  then  performing  the  integration 
over  tC     as  t  ->  oo  .   Although  the  surface  wave  decays  expon- 
entially in  time,  it  can  dominate  the  transient  term  in  some 
finite  time  interval.   The  parameter  which  determines  the 
extent  of  this  time  interval  is  the  effective  length  L  of  the 
region  in  which  the  plasma  equilibrium  is  inhomogeneous .   For 
the  profile  shown  in  Fig.  2,  L  =  (b  -  a)  .   As  L  ■*■  0,  the  plasma 
profile  approaches  one  having  a  sharp  boundary.   In  this  limit 
one  can  show  that  R  =  0  in  the  forced  component  of  v  while 
in  the  transient,  a,  y  and  u  .  are  identically  zero.   There- 
fore if  there  is  a  sharp  boundary,  the  transient  term  consists 
of  an  undamped  surface  wave  and  is  of  the  same  order  as  the 
forced  term.   This  behavior  can  be  explained  in  terms  of  the 
condensation  of  the  continuous  spectrum  of  the  inhomogeneous 
plasma  into  a  point  spectrum  as  L  -*■   0 . 
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The  transient  term  of  the  other  fields  which  characterize 
the  system  can  exhibit  a  time  behavior  which  differs  from  that 
of  v  ..      For  example,  by  differentiating  eq.  (88)  to  form  v  . , 
the  transient  of  v  ,  according  to  eq.  (25),  one  finds  that 
as  t  ■*■  oo  , 


dw 
Vzt  *  6(x)  ~dT  Cx)  sin  uaCx^  »  (89) 


i.e. 

Alfven  frequency.   The  transient  of  the  current  density  grows 

linearly  in  time, 


J  t  ^  Hx)E  -^  Cx)]n  cos  ua(x)t  .  C90) 


The  expressions  that  are  given  above  are  not  valid  at  the 

2     2 
point  x  where  co   =  oj  (x  )  .   The  reason  for  this  is  that  the 

singularities  of  Z   and  J   superimpose  to  form  a  singularity 

of  higher  order  in  y   .   Thus,  at  x  =  x  (u  )  , 

5m,  (w  -  co  ) 

v  (x  ,u)  =  Z  (x  ,oj)J  (to)  ^  2-   ,  (91) 

X   o*        vv  o'    S   '       10  -  OJ       ' 


which  implies, 


v    (x,t)   *   fa   t    cos   ooQt    ,  C92) 


as  t  ->  od  .   For  the  same  reason,  it  be  shown  that  at  x 


2k 


(93) 


Jy  *  t2  cos  u>Qt  ,  (.9*0 


The  temporal  growth  which  occurs  at  x  will  be  associated  In 
the  next  section  with  local  energy  absorption  at  that  point. 
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4.   Energy  Absorption 

In  this  section  we  shall  be  concerned  with  the  behavior 
of  the  plasma  energy  W  In  the  limit  t  -»■  oo  .   In  particular  we 
wish  to  determine  the  conditions  under  which  dW/dt  has  a  non- 
zero average  value  when  the  surface  current  varies  sinusoidally 
in  time.   If  the  average  of  dW/dt  over  the  period  of  the  applied 
frequency  is  positive,  then  one  can  conclude  that  the  plasma 
absorbs  energy. 

Equation  (33)  expresses  dW/dt  in  terms  of  v   and  p*  at  the 
plasma-vacuum  boundary.   In  the  limit  t  ■*■  <x>  ,  v  (b,t)  is  ob- 
tained from  eq.  (83), 


vx(x,t)  =  JoRv(b,wo)  sin  u>Qt  +jx,0>,w0)  cos  wQt,       (95) 


where 


iu^y    iw.  k>    iw~y   -k(c-b) 
V&,-0J  ♦  «£...„]  -  -  ^  e-«-b>  ♦  55a-  ^^   ,  (96) 

v  v         o 


while  p*(b,t)  can  be  derived  by  taking  the  inverse  Laplace 
transform  of  eq.  (43) } 

p*(b,t)    =   J   R    (b,w    )    sin  w  t    +   J  X  Cb,o)    )  cos   w   t  (97) 

^  opo  o  opo  o 


where 


Rp(bJ%)    +    iXp(b5.o)    =    By   el+lq]%)       •  (98) 
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By    combining  eqs.    (33),    (95)    and    (97),    the    following  expression 
for  the   time   rate   of   change   of   the   plasma  energy   is   obtained, 

(99) 

-2k(c-h) 
e 


dW    _         ^tL      J2    h»0yq1(M0)  _2k(c-b) 


l+kq(u    ) 
2 

k        2     ^vw»u,o/ilpv"»u,ov    "  "vvu>uVV 


2ttL. 


-|-   [Rp(b,0)o)Xv(b,u)o)    +   Rv(b,Wo)Xp(b,a)0)]sin    2wot    , 


where  q.(w  )  is  the  imaginary  part  of  q(u>  ).   According  to 

eq.  (99),  dW/dt  is  given  by  the  sum  of  a  term  which  is  constant 

in  time  and  terms  which  oscillate  sinusoidally  In  time  at  twice 

the  applied  frequency.   The  constant  term  Is  the  average  of 

dW/dt  over  the  period  T   =  27r/u   , 
o       o 

(|W\  =  ,^    *oWlU0)  je-2k(c-b)  (10Q) 

\dt/       k   |l+kq(a)0)|2  2 

and  the  time  varying  terms  are  due  to  the  energy  of  the  plasma 
fields  which  have  saturated  in  magnitude.  A  necessary  condition 
which  must  be  satisfied  in  order  that  q.(w  )  be  different  from 

zero  is  that  oo   lies  in  the  continuum,  i.e.,  there  must  exist 

2     2 
a  real  point  x  (u  )  in  the  plasma  at  which  u>  =  w  (x  ).  If  this 

condition  is  not  satisfied,  then  the  initial  conditions 

ijj  (0,u  )  =  0  and  4/(0, to  )  =  k  imply  ij>  is  real  for  x>  0  which  in 

turn  implies  that  q  is  real.   If  u   is  in  the  continuum,  then 

q.(w  )  can  be  expressed  in  terms  of  the  Wronskian  of  the 

independent  solutions,  f   and  f 2 ,  of  eq.  (40),  where  about 

x  (o>  )  ,  f -,  and  fp  are  given  by  eqs.  (60)  and  (6l),  and  the 
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Wronskian   is, 


•\,  df2         dri 


p(xo)g2CxoJu)0)    dco2 


(101) 


A(x,u)o)  dx      v   o 


CxJ       . 


Let  us  assume  that  <*>  (x)  is  a  monotonic  function  of  x  for 

0 _<  x  <   b,  so  that  there  exists  only  one  point  at  which 

a)   =  oj „(x).   Then  the  dependence  of  q.(w  )  on  W  is  established 

by  analytically  continuing  eq.  (62)  to  the  region  x>  x  (w  ) 

according  to  the  following  prescription, 


(102) 


xo(03o)    -  x    =   [x  -  xoCWo)]e 


-i€TT 


Kx,u>Q)    =    Ep(coQ)    -    i€Tr]g(x,a)o) 


+    [iU|xo(ojo)    -   x|g(x,u>0)    +    Cx-xo(a)Q)  )h(x,u)o)] 


where  €    =   sgn    (da)  /dx)      *      s,    and   constructing  iKx,(d    )/> '  (x,u    ) 
a         jc    (ui    j  o  o 

O        O 

This   yields, 


Im   C-jt)   =   -  €  tt- 


ftCx,w .  .) 


k*(x,o)0)| 


(103) 


where   use  has   been  made   of  eq.    (102)    with  9    =   -  tt.      Substitu- 
tion of   this   expression   in   eq.    (5*0    yields, 


k2B2  p(x    )g2(x„,a)n) 


yA^(b,w    )    U'(b,co    )r 


du 
d#<xo> 


(104) 
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2  ? 

where  use  was  made  of  the  identity  e  (doo  /dx)  =   du>  /dx  . 

a         a 

Since  qi  <  0,  eq.  (100)  implies  that  (dW/dt)   >  0  and  consequen- 
tly that  the  plasma  is  an  absorptive  medium  in  a  frequency  range 

which  coincides  with  the  continuous  spectrum. 

If  there  exists  more  than  one  real  x  (u>  )  at  which 

o   o 
2     2 
ooq  =  <d  (x),  the  analytic  continuation  defined  by  eq.  (102)  can 

be  performed  at  each  such  point  x  (x  ),  v  =  1,  2,  ...,  to 

obtain  q.  .   Specifically,  about  xV  ,  we  have, 


*  E  *v  ~=    Qv(a,o)[pv(a)o)flv(x'CV  +  f2v(x'Wo)]  '      (105) 
where , 

flv  =  gv(x'uo)  •  (106) 


f2v  =  ?vCx'Uo)  *"  Cxo~x)  +  (x"Xo)hv(x>%)   •         (107) 


g   and  b   are  analytic  at  xQ  and  normalized  such  that  they  are 

real  for  w   real.   Let  a  be  a  point  which  is  less  than  x  but 
o  v      ^  v 

in  the  domain  where  eq.  (105)  is  valid.   Then,  p   is  determined 
by  imposing  continuity  of  ip/i^1  at  x  =  a   ,  yielding 


where 


,  ,      f2v(av»M)  -  Vl  f2v<av»M) 
^  =  "  flv(av'w)  -  Vl  *Vav»w> 


(109) 
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To  obtain  Q  ,  we  impose  continuity  of  ij/  at  x  =  a   ,  but  since 
q(u>)  is  independent  of  Q   ,  its  form  will  not  be  of  concern 
to  us  here.   Analagous  to  what  appears  in  eq.  (102),  the 
analytic  continuation  of  i|i   from  x  <  xv  to  x  >  x  proceeds 
as  follows, 

(110) 
i|>v  =  Qv(u)[pv(«*))gvCx,u)  +  gv(x,03)^(xo-x)  +  (x-xQ)hv(x,w)], 

V 


Q  (u){[p  U)-l€  ir]g  (x,u)+g  U,u)ir.'(x-xJ)+  (x-x^)hv(x,u>)} 

x  >  xv 
o 

Qv(w){[pv(u))-  iev7r]?lv(x,w)  +^2v(x,a))}  ,   x  >  x^ 


2  ■«        ^ 

where  €   =  sgn  (do)  /dx)   v,  ^  ,    and  f,   and  f   are  analytic 

continuations  of  g  (x,w)  and  g  (x,u>)  S^  ^x-x0^  +  ^x_xo^  hv  ^XjU^ 

Let  p  .  be  the  imaginary 

part  of  p   .   Then  by  making  use  of  eq.  (110),  it  can  be  shown 


that, 


In) 


<k,  (x,<u) 


^  Cx  ,  u  ) 


tpvlo.)- V] a[^(x-:)-^(x'M)] — ,  ,(iiD 


I  [pv  (cjo)-iev7r]f^v  (x,u))  +  f^v  (x,a))  | 
where  from  eq.  (10  8)  , 


a[flv(av,M),  f2v(av,M)] 


(112) 
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a(v-l)i 


.  is  the  imaginary  part  of  a 


(v-1)  • 


The  Wronskians  of 


f     f   and  f\   ,  ^o  "have  the  same  functional  dependence 

lv      2v        -LV      dV 

on  x,  i.e., 


^(flv'  f2v}/ 

A(x,u)      dx  ^  o; 
Wlv.*2v>" 


(113) 


By    combining  eqs.    (109),    (111),    (112)    and    (113),    the    following 


recursion   relation   for   a   i   is    found 


where , 


0. 


€     T7$     p(x      ) 
V         VV        O 


;2(xV,cd) 
3v  v    o'     y 


dU)a        V 


-    a(V-l)i®v*Vp2(X0' 


dw 

—2-  fxul 


V  |A|l*l»^v-l*ivl 


<!<*>> 


v    ■    |A|[(p  .-It,,)*!.  *YSv  | 


1     i  2 


,    (114) 
(115) 

:=av+l 


Let  xn(w  )  be  the  zero  of  A  which  is  the  closest  to  the  plasma- 
vacuum  interface  at  x  =  b .   Then,  qi(w0)  is  obtained  by  sub- 
stituting eq.  (Ill)  into  eq.  (5*0  with  v  set  equal  to  n, 
yielding,  (116) 

,  (O    =    -it 


yA    (b,uQ) 


p(xg)sg(xg,o»0)  1^(^)1 

2 


l*n(b,»0)| 


iio' 


— ^  (xn) 
dx     ^o; 


kV 


e 


n"(n-Di  uA2(b>u   ) 


P2<*>g(y*'o>lQn(%>l 
l^tb,^)!2 


©. 


du 


^(xn) 
X       o 
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a,   -,>...  follows  from  the  recursion  relation  eq.  (114)  with 

Eq.  (114)  can  also  be  used 
to  show  that  each  singular  point  xv(u)  )  adds  a  negative 
definite  term  to  eq.  (116).   Therefore,  q. /   v  is  negative 
definite  which,  according  to  eq.  (100),  implies  uiW/dt^   >  0. 

Asymptotically  in  time,  the  energy  absorbed  by  the  plasma 
is  transferred  to  an  infinitesimally  narrow  region  about  x  (w  ) 
and  is  partitioned  so  as  to  increase  predominately  the  plasma 
motion  along  the  z  axis  in  this  region.   To  show  this  we  shall 
obtain  below  ^dW/dy  by  carrying  out  the  integration  over  x 
in  eq.  (34). 

It  was  pointed  out  in  the  previous  section  that  when  the 
applied  force  varies  sinusoidally  with  constant  amplitude  as 
t  ■+■  oo  except  at  x  =  x  (w  )  where  the  fields  grow  in  time. 
Therefore,  the  only  region  which  can  give  rise  to  a  non-zero 
average  of  dW/dt  is  some  narrow,  region  about  x   .   Moreover 
since  the  singular  eigenfunctions  of  v  are  square  integrable 
over  x,  the  contribution  to  W(t)  from  v  and  E   approaches  a 
constant.   Therefore,  to  leading  order  as  t  -*■  oo ,  we  obtain  with 
eq.  (3^), 


/dw\ 

\dt/ 


(117) 
x  (w  )+Ax  \ 

°   9  \  2ttL 

dx[|p(x)v2(x,t)+  ^-  k2B2(x)^(Xjt)]\_FX 


v  and  E  are  found  by  combining  eqs .  (42),  (51)  and  (52)  and 
then  performing  an  inverse  Laplace  transform.  For  x  %  x  (co  ) : 
this  leads  to,  in  the  limit  t  ■*■  oo  , 
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(118) 


du>2  J      sinkco    (x)-o)   ]t 

du^  J      sinkto    (x)-oo    ]t 

,z   «      |^(Uo)|g(xo).o)^(xo)-|  ^^  °        cos   „ot 

where  C  is  related  to  C  by, 


C(w)  =  C(u)Jo(a))  (119) 


Combining  eqs .  (117)  and  (118)  and  integrating  over  x  yields 
for  large  t, 


(Sf> 


2ttL      tt   £(io    )      g    (x    ,u)    )  rdco2  ,    J 

k  .2  pUo-)     dx  UoJ        2        •         U<^Uj 

K  CO  ' 

O 


It  can  be  shown  that  eq.  (120)  is  in  agreement  with  the 
expression  obtained  by  combining  eqs.  (100)  and  (104). 


■-'■• 


5.   Coil  Impedance 

In  an  experimental  situation  where  the  plasma  resonance 

is  excited  by  means  of  the  current  in  a  coil  surrounding  the 

plasma,  a  measurable  quantity  which  determines  both  the  power 

absorbed  by  the  plasma  and  the  energy  stored  outside  of  the 

resonant  region  is  the  coil  impedance  Z(w).   Since  the  current 

in  the  coil  is  effectively  a  surface  current,  the  impedance  is 

given  by  the  ratio  -6E  (to) /J  (oj),  where  J  and  6E  are  respective- 
y      S  S 

ly  the  Laplace  transforms  of  the  surface  current  and  the  vacuum 

electric  field  parallel  to  J  .   For  the  planar  model  in  Fig.  1, 

we  have, 

6E  Cc.to)  (121) 

Z  (to)  = V— 

j-CbO 

±   a3U_  e~kcslnh   kc    _    1   WVL  sinh   kb  -  kq  (u )  cosh   kb   e~k(2c-b) 


T  e         sinn   ice   -   i  f  i+kqM 


where  6Er(c,w)  is  expressed  in  terms  of  J  (w)  by  eqs .  (44)  and 
y  o 

(53).   When  J  (t)  varies  sinusoidally  in  time  in  the  manner 
shown  by  eq.  (82),  the  average  power  abosrbed  by  the  plasma  as 
t  ^  cd  is  realted  to  Z  as  follows, 


(122) 


e-2k(c-b)  (123) 


Equations  (100)  and  (122)  are  seen  to  be  in  agreement.   The 
imaginary  part  of  Z  determines  the  total  average  energy  in  the 
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2ttL 

J2 
■  Re[Z(u)o)]    -f 

re    from  eq.     (121) , 

Re[Z(a))]    =    - 

coyq1(a)) 

|1+   kq(aj)  |2 

plasma  and  vacuum  regions  outside  the  resonance  zones.   To 
illustrate  in  what  manner  Z  typically  depends  on  w,  let  us 
consider  the  piece-wise  smooth  profile  in  Fig.  2.   For  this 
case,  q  has  an  explicit   representation  in  terms  of  Bessel 
functions  that  can  be  obtained  by  combining  eq.  (5^)  and  (67) 
where  0  is  to  be  set  equal  to  -it  in  eq.  (67).   This  yields, 


[p(co)-iTr]Io[k(b-xo)]-Ko[k(b-xo)] 
kqU)  =  [P(oo)-iTrJl1Lk(b-xo)]+K1Lk(b-xo)J   ■  (124) 


where  p(w)  is  given  by  eq.  (58)  and, 

?           2    2  2    2                   (125) 

up-kBVy  k   B/y 

x    -  a   =   -5 =-* — y? (b-a),      b-x      =   -* * — ~ — * (b-a) 

u>2pp+k2(B2-B2)/p  °        u>2pp+k2(B2-B2)/y 

It   can  be    shown  with  eqs .    (58)    and    (124)    that    q   is    complex   in 

2    2 
the   continuum  w   >   to    (=    k  B  /p   u ) ,      but   real    for   0  <  to  <  w   . 

The   imaginary   part   of  q  which   is   to  be    substituted   into   eq. 

(123)    is    given  by, 

(126) 

,  ,      \     _  TT  _    _1 

qi   w  Hb-XQ)     |[p(U3)_i7T]ii[k(b_Xo)]+K1[k(b-xo)]|2 

The  real  and  imaginary  parts  of  Z  are  plotted  in  Fig.  6 
as  a  function  of  f(=  w/2tt)  for  a  hydrogen  and  a  deuterium 
plasma.   The  parameter  values  chosen  are  typical  of  a  high 
beta,  high  density  plasma.   The  real  part  of  Z,  Zr,  is  identi- 
cally zero  up  to  the  low  frequency  end  of  the  continuum  at 
f  =  f  .   For  f  >  f   ,  Z   peaks  at  about  2  MHz  for  a  deuterium 
plasma  and  at  2  •  8  MHz  for  a  hydrogen  plasma,  and  then 


y^ 


approaches  zero  as  f  ■*  oo  .   In  the  vincinity  of  the  peak,  the 
form  of  Z(w)  suggests  the  existence  of  a  damped  resonance. 
This  resonant  behavior  is  due  to  the  presence  of  a  complex  zero 
of  the  denominator,  1  +  kq(oj),  of  eq.  (121)  near  the  real  u  axis. 
This  zero  is  the  complex  frequency  of  the  damped  surface  wave 
which  was  discussed  in  section  3. 

Fig.  6  shows  that  the  resistive  component  of  the  coil 
impedance  has  a  peak  value  of  about  one  Ohm.   To  estimate  the 
efficiency  of  the  absorption  process,  one  should  compare  this 
value  with  the  surface  resistance  of  the  coil.   For  example, 
the  surface  resistance  of  copper  wire  at  2  MHz  is  about 
0.0007   Ohm,  implying  that  if  the  coil  is  made  of  copper,  only 
0.07   per  cent  of  the  power  delivered  to  the  coil  would  be  lost 
to  Joule  heating  of  the  wire,  while  93  per  cent  of  the  power 
would  be  abosrbed  by  the  plasma. 
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Plasma  Configuration. 
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FIG.  2 


Piece-wise   smooth  equilibrium  profiles.      In  the   vacuum 
region    |x|    >  b,    a  sheet   current   is    located   at   x  =  ±  c. 
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Branch  points    (x )    and  branch  cuts    (aaa)    of  ^(x,ao) 
for   the   piece-wise    smooth  profile   in  Fig.    2. 
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FIG.  4 

Branch  points  (x )  and  branch  cuts  (<w\)  of  ^(x,oo)  (Fig.  (b)) 

2 
for  the  case  where  co  depends  on  x  in  the  form  shown  in 

Fig.  (a). 
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FIG. 5 


The  contour  £  along  which  the  inverse  Laplace  transform 
is  to  be  accomplished  and  the  contour  £'  obtained  by 
deforming  £  around  the  singularities  of  v  (x,od). 
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FIG.  6 


Coil  impedance  Z  vs.  frequency  for  a  hydrogen  (H  )  and 

a  deuterium  (D+)  plasma.   The  plasma  profile  is  shown  in 
Fig.  2  with  the  following  parameter  values  assumed: 

2   2 
a  =  2cm,  b=  2.4cm,  c=2.6cm,  B„  =  15kg,  B_/B„=0.5J 


15kg,  B^/B* 


k=  0.5cm        and  n  =  10      p/cirr    where    n   is    the    equilibrium 
particle   density. 


44 


References 

Baldwin,  D.E.  and  Ignat ,  D.W.  1969  Phys,  Fluids,  12,  697- 
Barston,  E.M.  196M  Ann.  Phys.29_,  282. 
Boris,  J. P.  1968  Ph.  D.  Dissertation,  Princeton  Univ. 
Brownell,  J.H.  1972  Physica,  6l_,  289. 

Goedbloed,  J. P.  and  Sakanaka,  P.H.  197^  to  be  published  in 
Phys.  Fluids. 

Grad,  H.  1969  Phys.  Today,  22,  3H. 

Grad,  H.  1973  Proc.  Nat.  Acad.  Sci.  USA  70_,  3277- 

rd 

Grossmann,  W.,  Kaufmann,  M.  and  Neuhauser  J.  1973  3 —  Int. 

Symp.  of  Toroidal  Plasma  Confinement,  Garching,  paper  E3- 

Grossmann,  W.  and  Tataronis,  J. A.  1973  Z.  Physik,  26l_,  217. 
Hasegawa,  A.  1973  Private  Communication. 
Landau,  L.  19^2  J.  Phys.  USSR  10,  25. 

Puri,  S.  1973  Report  IPP  IV/59,  Max-Planck-Institut  fur 
Plasmaphysik,  Garching. 

Sturrock,  P. A.  1961  M.W.  Lab.  Report  784,  Stanford  Univ., 
Stanford,  Calif. 

Tataronis,  J. A.  and  Grossmann,  W.  1972  Proc.  of  The  Second 

Topical  Conference  on  Pulsed  High-Beta  Plasmas,  Garching, 
IPP-Report  1/127,  papers  B-5  and  B-6. 

Tataronis,  J. A.  and  Grossmann,  W.  1973  Z.  Physik,  26l,  203- 


45 


Uberoi,  C.  1972  Phys .  Fluids,  15,  1673- 

Watson,  G.N.  1958  A  Treatise  on  the  Theory  of  Bessel  Functions 

(2nd  ed.),  Cambridge  Univ.  Press. 
Weitzner,  H.  1963  Phys.  Fluids,  6,  1123. 
Weitzner,  H.  1973  Phys.  Fluids,  16,  237. 


46 


This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  Implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  In  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
In  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  Information,  apparatus, 
method,  or  process  disclosed  in  this 
report. 

As  used  in  the  above,  "person  acting  on  behalf 
of  the  Commission"  includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 


47 


MAY  ?0  1<?74 

DATE  DUE 

m-r. 

-. 

I 

™,«»™«.... 

COO-3077-;46 
Tataronis 


c.l 


COO-3077-46  c.l 

COO-3077-46  c#1 

Tataronis 


JjJLenergy   absoretiondue   to 


the    cnnt.lminip    spectnim. 

BORROWER  S    NAME 


^l^^rfe^^^ 


N.Y.U.  Courant  Institute  of 
Mathematical  Sciences 

251  Mercer  St. 
New  York,  N.  Y.    10012 


